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1.  Introduction. 


System  theory  is  concerned  with  the  modelling  and  analysis 
of  phenomena  both  natural  and  man-made.  It  is  a discipline  whose 
formal  beginnings  go  back  at  least  to  Watt  and  Maxwell  and  much 
of  its  motivation  stems  from  engineering  problems.  Before  World 
War  II,  a system  design  and  analysis  were  primarily  an  art.  During 
and  after  the  war,  techniques  based  on  complex  variable  theory  were 
developed  and  applied  primarily  to  single  input,  single  output 
systems  represented  by  a rational  function,  called  the  transfer 
function.  The  theory  of  servomechanisms  developed  rapidly  from 
the  end  of  the  war  to  the  early  fifties  and  time-domain  methods  were 
applied.  The  representation  of  transfer  functions  via  linear, 
constant-coefficient,  differential  equations  led  to  a renewed 
interest  in  so-called  state  space  methods.  The  rapid  development 
of  the  theory  followed  and  continues  today.  However,  the  increasing 
complexity  of  the  engineering  and  economic  problems  considered 

l 

required  greater  mathematical  sophistication.  Little  did  one  of 
the  authors  dream  while  he  was  a student  of  Professor  Zariski,  that 
the  techniques  of  modern  algebraic  geometry  could  (and  would)  be 
applied  in  a critical  and  essential  way  in  system  theory.  Our 
purpose,  in  this  paper,  is  to  illustrate  several  of  the  many  applica- 
tions of  algebraic  geometry  to  linear  system  theory.  These 
applications  are  both  mathematically  and  practically  non-trivial. 

We  begin,  in  Section  2,  with  an  analysis  of  two  (classical) 
questions;  namely,  (i)  when  is  a meromorphic  function  of  the  form 

oo 

_ _ £ 

f(s)  = L h.s  rational?,  and,  (ii)  when  is  such  a rational  function 
£ = 1 * 


* 


4 

f 


j 


stable  (i.e.  has  poles  in  the  left-half  plane)?  The  first  question 

was  answered  by  Hankel  ([20])  and  the  second  was  answered  by  Cauchy, 

Hermite,  Routh  and  Hurwitz  ([15],  [29],  [30]).  Our  treatment  serves 

to  motivate  the  subsequent  development.  The  concept  of  a time- 

invariant,  finite-dimensional  linear  system  is  introduced  in 

Section  3.  Such  systems  admit  an  external  description  as  either  a 

matrix,  T(s),  of  rational  functions  or  as  an  input-output  map 
00 

f(s)  = l L0s  with  the  L0  constant  matrices  and  rank  < °° 

£=1  * * 1 

where  = (L^  + j ^)  is  the  Hankel  matrix  of  f.  Corresponding  to 
each  external  description  is  an  internal  description  or  representation. 
Thus,  a pair  of  polynomial  matrices  (R(s),P(s))  realizes  T(s) 
if  T(s)  = R(s)P  *(s)  and  R(s),P(s)  are  relatively  right  prime. 

The  group  <2r  of  polynomial  matrices  with  determinant  a unit 
acts  on  such  realizations.  For  an  input-output  map  f,  a triple 
(A,B,C)  of  constant  matrices  represents  f if  Hj  = (CA1  + -I’1B). 

The  general  linear  group  G (of  appropriate  dimension)  acts  on 
such  triples  via  g-(A,B,C)  = (gAg’1 .gBjCg’1) . Three  main  topics 
are  treated  in  Section  4.  The  first  is  the  construction  of  a 
moduli  space  for  linear  systems  under  the  action  of  G and  a proof 
of  a theorem  on  "realization  with  parameters",  based  on  Zariski's 
Main  Theorem  and  Mumford's  geometric  invariant  theory.  It  turns 
out  that  the  "exceptional  points"  are  precisely  the  representations 
that  are  undersirable  from  a practical  viewpoint.  Thus,  the  mathe- 
matical and  physical  considerations  are  exactly  congruent.  The 
second  ;opic  is  a brief  study  of  the  biregular  and  topological 
invariants  of  the  moduli  space  with  an  emphasis  on  the  impact  these 
invariants  have  on  system  theoretic  questions  analogous  to  the 


questions  raised  in  Section  2.  The  third  topic  is  the  considera- 
tion of  moduli  spaces  for  systems  with  richer  symmetries  such  as 
Hamiltonian  or  symmetric  systems.  In  fact,  we  sketch  the 
equivalence  between  an  old  problem  in  network  theory  ([53])  and 
what  might  be  called  the  "serre  conjecture  for  quadratic  modules". 
Recent  results  on  quadratic  modules  ([1])  are  precisely  what  is 
needed  for  the  solution  of  the  network  theory  problem.  The  vital 
concept  of  feedback  is  considered  in  Section  5 together  with  its 
relation  to  the  problem  of  pole  placement  (stabilization)  or  co- 
efficient assignment.  Equivalence  under  feedback  involves  the 
action  of  a group  which  is  not  reductive.  Nonetheless,  a moduli 
space  can  be  constructed  and  its  invariants  calculated.  For  con- 
trollable systems,  the  result  can  be  found  in  [7],  [28],  [32]  and 
[42] . The  general  case  is  completely  analyzed  in  [19]  . It 
should  also  be  possible  to  extend  Mumford's  geometric  invariant 
theory  to  groups  of  feedback  type.  A critical  portion  of  the 
feedback  invariants  is  the  Kronecker  set  of  indices  (k^,...,k  ) 
where  the  are  non-negative  integers.  This  set  was  used  to 

define  an  ordering  for  systems  by  Rosenbrock  ([42]).  It  turns 
out  that  this  ordering  is  precisely  the  Harder-Narasimhan  ordering 
C[45])  and  we  conclude  with  an  indication  of  the  use  of  this 
complex  of  ideas  in  system  theoretic  problems.  We  hope  to  indicate 
throughout  the  richness  of  the  relationship  between  algebraic 
geometry  and  system  theory  both  in  terms  of  the  application  of 
algebraic  geometry  to  system  theory  and  in  terms  of  the  generation 
of  problems  in  algebraic  geometry  from  the  practical  considerations 
of  system  theory. 


2 . Routh-Hurwitz  Theory 


We  begin  by  considering  two  questions,  which  were  posed  and 
solved  in  the  19^  century  by  Cauchy,  Hankel,  Hermite,  Hurwitz, 
Maxwell,  Routh  and  others,  and  which  now  form  part  of  the  folklore 
of  linear  system  theory. 

Question  1.  When  is  a proper  (i.e.,  vanishing  at  00 ) meromorphic 
function  f on  C rational? 

Question  2.  If  such  an  f is  real,  when  do  its  poles  lie  in  the 
left-half  plane? 

In  fact,  we  can  ask  the  same  questions  for  proper,  matrix- 
valued meromorphic  functions  and,  as  a motivating  example,  consider 
the  resolvent  or  transfer  function,  f(s)  = (sI-F)  \ of  a differential 
equation,  (^-  - F)x  = 0.  This  is  still  more  interesting  in  the 
control  system  case  (see  Section  3),  but  serves  to  explain  Maxwell's 
interest  ([38])  in  such  questions.  Question  1 can  be  thought  of 
in  terms  of  an  algebro-geometric  classification  of  differential 
equations,  whereas  Question  2 characterizes  asymptotic  stability. 

One  approach  to  Question  1,  due  to  Hankel  ([20]),  is  to  form 
the  matrix,  = (h^+j_^),  where  h^  is  given  by 

CD 

f(S)  = l h s'*.  (2.1) 

£=1  * 

Hankel1 s Theorem  then  asserts:  f is  rational  if,  and  only  if, 
rank  * n < 00 , as  one  can  plainly  see.  In  this  case,  f induces 
a holomorphic  map 


2-2 


f: 


(2.2) 


whose  degree  is  given  by 


deg^f  = rankfH^.  (2.3) 

Now,  in  case  £ is  real,  we  can  restrict 

f:  1R1P  1 -*■  1R1P1  (2.4) 

and,  since  £ K Hf  is  bijective,  we  can  ask  for  deg^  £ in  terms 
of  H£.  In  this  way,  we  obtain  the  beautiful  result  which  underlies 
the  Routh-Hurwitz  conditions: 

Theorem  2.5  (Hermite-Hurwitz)  deg^  f = signature(H£) . 

This  was  discovered  by  Hermite  in  1856  ([29  ])  in  case  £ has 
distinct  poles,  and  later  extended  by  Hurwitz.  Of  course,  deg  ^f 
was  expressed  in  a different  way: 

Definition  2,6  (Cauchy)  The  local  index  of  a real,  rational  f at 
a real  pole  Xq  is  +1  if  f changes  from  +00  to  -00 , -1  if 
the  alternate  change  occurs,  and  0 if  f does  not  change  sign 
while  passing  though  Xq.  The  index  of  f,  C(f),  is  the  sum  of  the 
local  indices. 

C(f ) , which  is  clearly  the  winding  number  of  the  map  in  (2.4), 
was  defined  by  Cauchy  in  [ 15] . In  part  I of  [ 15] , he  uses  the 
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Cauchy  index  to  compute  the  number  of  real  roots  of  a real  poly- 
nomial (generalizing,  among  other  things,  Descartes'  rule  of  signs), 
the  number  of  negative  real  roots,  and  related  questions.  In  part 
II,  he  uses  the  Cauchy  index  to  define  and  evaluate  the  local  index 
of  a non-degenerate  plane  vector  field  at  an  equilibrium  point.  Of 
course,  this  latter  application  was  generalized  by  Kronecker,  for 
n > 2,  and  others.  On  the  other  hand,  the  computation  of  the  number 
of  negative  real  zeroes  of  a polynomial  is  a bit  harder  than 
Question  2 and  this  computation  was,  in  fact,  the  starting  point 
for  Hurwitz  in  [30].  Following  [30],  we  assume  for  the  moment 

that  f has  no  poles  on  the  imaginary  axis.  Thus,  if  g is  the 
monic  denominator  of  f,  having  no  roots  in  common  with  the 

associated  numerator  h,  computing  the  change  in  arg  g(-is),  s between 
+“  and  -®»,  is  equivalent  to  computing  p - q,  the  number  of  zeroes 
in  the  left-half  plane  minus  the  number  of  zeroes  in  the  right- 
half  plane.  This  is,  by  trigonometry,  the  Cauchy  index  of  the 
rational  function  v/u  where 

cg(-is)  = (u+iv) (s) , 

and  c is  a complex  constant  rendering  v/u  proper.  By  the 
Hermite -Hurwitz  Theorem, 

p - q = sign  Hu/y.  C2.7) 


Thus,  by  the  Jacobi -Frobenius  algorithm,  one  obtains  numeric 
criteria  (in  terms  of  the  coefficients  of  g!)  for  Question  2;  i.e., 
for  g to  be  a Hurwitz  polynomial.  Explicitly,  this  is 
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positive-definiteness  of  a quadratic  form,  and  one  should  stress 
that  it  is  a known  fact  in  the  folklore  that  Hu/V  is  intimately 
related  to  a Liapunov  function  for  = 0.  In  fact,  this  is 

very  much  the  spirit  of  Hermite's  treatment. 

We  close  the  section  by  sketching  a moduli- theoretic  proof  of 
Theorem  2.5,  due  to  R.  W.  Brockett  ([  2 1).  With  the  notation  as 
above,  consider 

Rat(n)  = {(g,h):  g/h  is  proper  and  Cg,h)  = 1 in  IR[s]}. 

As  the  colocus  of  the  resultant,  Rat(n)  is  Zariski  open  in  IR^n, 
the  Euclidean  topology  being  the  same  as  the  compact-open  topology 
on  the  maps  (2.4).  In  particular,  the  Cauchy  index  p - q is  constant 
on  Euclidean  components  of  Rat(n)  and  there  are  at  least  n + 1 
of  these.  By  continuity,  signCH^)  is  constant  on  components  of 
Rat(n),  in  light  of  (2.3).  Thus,  already,  a general  position 
argument  reduces  Theorem  2.5  to  the  result  obtained  by  Hermite. 

Indeed,  if  Rat(p,q)  is  the  submanifold  on  which  the  Cauchy  index 
is  p - q,  then,  by  a somewhat  tedious  but  intuitive  canonical 
form  argument  based  on  divisors,  Brockett  proved 

Theorem  2.8.  Rat(p,q)  is  connected. 

It  is  entirely  trivial  to  evaluate  the  induced  bij ection  between 
values  of  the  Cauchy  index  and  the  signature  of  H^.  This  yields 
Theorem  2.5  as  well  as  some  insight  into  the  global  structure  of 
Rat(n)  (i.e.,  of  scalar  input-output  linear  systems).  It  is  very 
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natural  then  to  study  the  topology  of  the  submanifolds,  Rat(p,q)c 
fip -q(S^).  For  example,  based  on  unpublished  work  of  R.  Brockett 
and  G.  Segal,  one  can  show  that  each  Rat(p,q)  can  be  given  the 
structure  of  a Stein  manifold  and,  in  articular,  deduce  the 
vanishing  of  its  higher  cohomology  groups.  Their  work  also  shows 
that  Rat(p,q)  is  simply  connected  only  in  case  either  p = n, 
or  q = n,  in  which  case  Rat(n,0)  = Rat(0,n)  ~ ]R2n.  It  is  also 
known  that  Rat(n-l,l)  = Rat(l,n-1)  * S1  xR211'1  (see  [2]).  In 
another  direction  ([13]),  from  all  of  the  above,  we  can  construct 
a vector  bundle  V+  on  Rat(p,q),  obtained  by  assigning  the 
positive  eigenspace  of  to  f e Rat(p,q).  This  is,  in  general, 

non-trivial.  For  example,  on  Rat(l,l),  V+  induces  a classifying  map 

T+  : Rat  (1 , 1)  - 1RIPL 

which  is  not  homotopic  to  a constant  map.  Indeed  ([2]),  under 
the  isomorphism  Rat(l,l)  = S^  * ]r\  T+  is  just  projection  on 
the  first  factor.  In  section  four,  we  will  interpret  these 
obstructions  in  a system- theoretic  context. 


3 . Linear  Systems 


The  class  of  systems  which  we  consider  is  the  class  of  finite- 
dimensional, linear  time- invariant  multivariable  systems.  Such 
systems  can  be  characterized  via  an  external  intrinsic  description 
using  either  the  transfer  matrix  or  an  input-output  mapping. 

Associated  with  each  description  is  an  internal  representation  and 
an  appropriate  notion  of  equivalence  which  is  defined  by  the  action 
of  an  algebraic  group.  The  critical  link  between  these  descriptions 
is  provided  by  ( a generalization  of)  Hankel's  Theorem  on  the 
rationality  of  proper  matrix  valued  "meromorphic"  functions. 

Now  let  R be  a Noetherian  integral  domain  which  is  integrally 
closed  in  its  quotient  field  K and  let  R[s]  be  a polynomial  ring 
in  s over  R.  We  call  an  element  n(s)/d(s)  of  K(s)  proper 
if  degree  n(s)  < degree  d(s) . Let  Mp  m(0  denote  the  set  of 
p * m matrices  with  entries  in  • . Then  we  have: 

Definition  3.1.  Let  Ep  m = {T(s)  £ Mp  m(K(s))|  the  entries  t^  (s) 
of  T(s)  are  proper).  Elements  of  Ep  m are  called  proper  transfer 
matrices . 

Definition  3.2.  Let  Tfs)  be  an  element  of  E _ . A pair 
p,m 

(R(s),P(s))  with  R(s)  e M (Kts))  and  P(s)  an  invertible  element 

p , m 

of  Mm  m(K[s])  is  called  a representation  of  TCs)  if 

T(s)  = R(s)P  X(s) . A representation  is  called  a realization  (or 

minimal  representation)  of  T(s)  if  R(s)  and  P(s)  are  relatively 


right  prime. 

If  T(s)  e E . then,  letcing 

p ,m 

A (s) 

be  the  (monic) 

least 

common  multiple  of  the  denominators 

of  the 

entries  t^  (s) 

of  TCs), 
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we  have  T(s)  = R(s)P  1 C s ) where  R(s)  = TQ  + T^s  + •••  + T^^s11 

P(s)  = A(s)Im  and  so,  representations  exist.  Moreover,  let 

= (U(s)  e M m(K[s]):  det  U(s)  is  a unit  i.e.,  a non-zero  ele- 
m m,m  — 

ment  of  K }.  Note  that  <2^  is  an  algebraic  group  which  acts  on 
pairs  (R(s),P(s))  via  right  multiplication  and  preserves  the 
property  of  being  relatively  right  prime.  We  now  have: 


Theorem  3.3.  Given  Tfs)  £ E , there  exist  realizations  of  T(s) 

v p,m 

and  any  two  realizations  are  equivalent  under  the  action  of  <2^. 

(For  a proof,  see  [18],  [49].) 

Let  T(s)  be  an  element  of  E and  let  o~(s)  be  an  element 

p,m  I 

of  Mp+jn  m(K[s])  which  corresponds  to  T(s).  In  other  words,  Oj,(s) 

]~RT(sf 

is  an  element  of  M m(K[s])  such  that  cr _(s)  = with 

P m,m  “ T P_(S) 


R^.(s),P^.(s)  relatively  right  prime  and  T(s)  = Rt(s)Ptls).  Any 

such  oT(s)  is  called  a linear  system  with  transfer  matrix  T(s). 

If  E is  viewed  in  this  way,  then  E is  stable  under  the 

p,m  p,m 

action  of  ^ and  the  orbits  correspond  to  the  transfer  matrices, 
m 

More  precisely,  if  S „ CM  x fK[s])  is  the  set  of  all  linear 
r ' p»m  p+m,m  — 

systems,  then  S m is  stable  under  the  action  of  and  the 

' 7 p,m  m 

transfer  matrix  is  a complete  invariant  for  this  action.  These 
results  can  be  interpreted  in  still  another  way  which  is  frequently 
useful  in  applications.  Namely,  let  cr(s)  be  an  element  of  m 

and  let  ck  (s)  be  the  j t^1  column  of  o(s).  Then  o^Cs) am(s) 

are  free  over  K[s]  and  we  let  Mg  be  the  free  module  with 
generators  cf^  (s)  , . . . ,cr  (s)  . MQ  is  called  the  system  module  of 
o(s)  and  is  also  a complete  invariant  for  the  action  of  % . 
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, then  na  = degree  det  Pa(s)  is  called  the 

[ 

McMillan  degree  of  o(s)  and  nQ  is  also  an  invariant  for  the 

action  of  °ii  . It  is  rather  tempting  to  heuristically  view 
m 

(R0(s) ,P0(s) ) as  "homogeneous  coordinates"  relative  to  multiplica- 
tion by  U(s)  in  <2'm.  Indeed,  in  the  scalar  case  (p  = m = 1) , 
this  is  expressed  in  (2.2).  More  precisely,  cr(s ) induces  a map  T0 
of  IP1  into  GrassK(m,m+p)  given  by 

T0(s)  = ((Ra(s)w,Pa(s)w|  w e Km>  (3.4) 

for  s f 00  and  by  Km  for  s = °°.  is  essentially  the  graph 

of  T0(s)  by  virtue  of  the  fact  that  Ra(s) ,Pa(s)  are  relatively 
right  prime  ( [27]  ) . In  addition,  if  K is  algebraically  closed, 
then  (2.3)  generalizes  also  to  the  following 

de^K~a  = na  (3'5) 

which  was  first  proved  over  C by  Hermann  and  Martin  ([27]).  We 

shall  soon  see  that  n is  also  the  rank  of  a Hankel  matrix. 

o 

We  now  turn  our  attention  to  the  alternate  external 
description  of  linear  systems.  Consider  the  (formal)  Laurent  series 

00 

f(s)  = l L.s'*  (3.6) 

1=1  % 

where  e M^  m(K).  We  associate  with  f the  (generalized) 

Hankel  matrix 


If  o(s) 


PQ(s) 


Hf  - 


(3.7) 


and  we  have: 


Definition  3.8.  f is  admissible  if  rank  = n^  < ®.  Let 

£ f 

E „ = (f(s)l  f is  admissible}.  Elements  of  E are  called 

p,m  1 p,m 

proper  input-output  maps  and  is  called  the  dimension  of  f. 


Definition  3.9.  Let  f(s)  be  an  element  of  Er  . A triple 
p,m  * 

(A,B,C)  with  A e M (K) , B e M fK)  , C e M (K)  is  called 
’ n,n  — n,m  — ' ’ p,n  — 

a representation  of  f(s)  if  CA  B = L^  for  l = 1 A 

representation  is  called  a realization  (or  minimal  representation) 
of  f(s)  if  n = n^. 


Definition  3.10.  A representation  (A,B,C)  of  f(s)  is  called 
controllable  if  rank  5f(A,B,C)  = rank  [B ,AB , . . . ,An’ ^B]  = n and 
a representation  (A,B,C)  of  f(s)  is  called  observable  if 
rank  <?(A,B,C)  = rank [C  \ A ' C ' , . . . , (A ' )n_1C ' ] = n. 


If  f(s)  is  a proper  input-output  map,  then  it  is  well-known 
that  representations  and  realizations  exist  ([33]).  Moreover,  if 
(A,B,C)  is  a representation  of  f(s),  then  C(sl-A)  *B  = T£(s) 
is  an  element  of  E^  m and  (A,B,C)  can  be  called  a state-space 
representation  of  T^s).  Similarly,  if  (A,B,C)  is  a realization 
of  f(s) , then  C(sl-A)  ^B  = Tf(s)  is  an  element  of  E and 
(A,B,C)  is  called  a state-space  realization  of  T^s).  Thus, 

f 

there  is  a natural  mapping  of  E^  m into  E^  m and  it  can  be 
shown  that  this  mapping  is  bijective  by  virtue  of  Hankel's 
Theorem  ([20]). 
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Now,  if  (A,B,C)  is  a triple  with  A £ M (K) , B e M (K) 

IL  y II  _ II  y III 

C e M (K)  , then  GL(n,K)  acts  on  (A,B,C)  in  the  following 
p , m 

way:  (A,B,C)  ->  (gAg  ^gB.Cg  , g e GL(n,K).  It  is  clear  that 

controllability  and  observability  are  preserved  under  this  action. 

We  now  have: 

Theorem  3.11.  Given  f(s)  e m (or,  equivalently,  T^(s)  e m) , 
there  exist  realizations  of  f(s)  (state-space  realizations  of 
Tf(s))  and  any  two  realizations  are  equivalent  under  the  action  of 
GL(n£,K) . Moreover,  any  realization  of  f(s)  is  both  controllable 
and  observable.  Finally,  if  (A,B,C)  and  (A^B^jC^)  are 
realizations  of  f(s),  then  the  g e GL(n£,K)  such  that  (gAg  1,gB,Cg  1) 
= (A^B^C^)  is  unique. 

(For  a proof,  see  [ 3 1 , [17]). 

The  final  part  of  Theorem  3.11  is  often  referred  to  as  the 
state-space  isomorphism  theorem. 

So,  let  f(s)  e and  let  (A,B,C)  be  an  element  of 

2 P 9 ™ 

nt+nf(m+p)  2 r , 

A (K)  = /An  n(-  p-'  which  corresponds  to  f(s).  In  other 

words,  (A,B,C)  is  an  element  of  An  +n(m+P)  such  that  (A,B,C) 

is  both  controllable  and  observable  and  T^s)  = C(sl-A)  and 

00 

f(s)  18  I (CA^'"1B)s"A.  Any  such  triple  (A,B,C)  is  called  a 
A*1 

linear  system  with  input-output  map  f(s).  Let  m be  the  set 

of  all  linear  systems  of  dimension  n.  Then  m is  stable  under 
7 p,m 

the  action  of  GL(n,K)  and  the  orbits  correspond  to  the  proper 

input-output  maps  of  dimension  n.  In  the  next  section,  we  shall 

characterize  the  moduli  space  of  Sn  „ under  this  action.  Finally, 

r p,m 

we  remark  that  the  dimension  of  f(s)  is  precisely  the  McMillan 


• »«■«>  HHM 


degree  of  T^(s) . 

So  far,  we  have  considered  notions  of  equivalence  based  on 
classical  groups  such  as  and  :&L(n,K).  However,  in 

Section  5,  we  will  examine  the  critical  concept  of  feedback 
equivalence.  This  notion  requires,  amongst  other  things,  an 
additional  set  of  invariants  known  as  the  Kronecker  indices  ([14], 
[32])  which  we  now  introduce  in  two  ways. 

If  P(s)  e Mm  m(K[s]),  then  let  Kj (P)  = max  degree (p ^ (s) 
i = l,...,m}  be  the  j column  degree  of  P(s)  so  that  P(s)  can 
be  written  in  the  form 


PCs)  = Ac(P)diag[s  1,...,sn]  + P-^s) 


(3.12) 


where  k.(P^)  < k. (P)  and  A (P)  e M (K) . P(s)  is  called 
j 1 j c m,m  — 


column  proper  if  Ac(P)  e GL(m,K).  If  o = a( s)  = 


Pa(s) 


is  an 


is  an  element  of  S , then  it  is  well-known  that  there  is  a 

P »m 

U(s)  e such  that  Pa(s)U(s)  is  column  proper  ([  ]).  We  now 

have : 


Definition  3.13.  Let  o be  an  element  of  S „ and  let  P(s)  be 

p,m 

a column  proper  element  of  m(K[s])  such  that  P(s)  = Pa(s)U(s) 
for  some  U(s)  e Then  the  set  of  integers  <a  = ( (P)  , . . . , Km(P)  ) 

is  called  the  Kronecker  set  of  o. 


Theorem  3.14.  Let  a be  an  element  of  S 


Then  (i)  is 


well-defined,  and  (ii)  if  t = aU  for  some  U e ^m,  then 
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(as  sets) . 


(i.e.,  <a  is  an  invariant  for  the  action  of 


Alternatively,  in  state-space  form,  let  (A,B,C)  be  con- 
trollable so  that  rank [B , AB , . . . , An  ^B]  = n.  Then  it  is  possible 
to  lexicographically  order  the  first  n linearly  independent 

columns  of  [B,AB,...,An  ^B]  as  follows:  b, ,Ab, , . . . ,A  b,  , 

<2-l  Km-i  11  -1 

b2,...,A  b2»...,A  bm  where  b^,...,bm  are  the  columns  of  B? 

The  set  of  integers  k = (k^,..,.,k  ) is  called  the  Kronecker  set  of 

(A,B,C).  It  can  be  shown  that  < is  an  invariant  under  the  action 

of  GL(n,K)  ([7  ]).  For  example,  if  we  consider  the  pencil 

[A-sI,B],  then  equivalence  under  GL(n,K)  induces  a strict 

equivalence  (in  the  classical  sense  of  Kronecker  [ ] ) of  pencils 

and  the  integers  (k^ , . . . , k ) are  a complete  invariant  for  strict 

equivalence  of  pencils.  Explicitly,  following  Kronecker,  the  < ^ 

may  be  computed  by  the  process:  choose  <p^(s)  e Ker  [A-sI  ,B]  ® K 

of  minimal  degree  in  s^  then  choose  e Ker  [A-sI  ,B]  ® K, 

independent  of  4>1(s),  and  of  minimal  degree  <2  among  all  such 

<Ks) , etc.  (cf.  (10  ]).  We  also  have: 


Theorem  3.15.  Let  (A,B,C)  be  a state-space  realization  of  T(s) 
and  let  cr^,(s)  be  a realization  of  T(s).  Then  the  Kronecker  sets 

are  the  same. 

(For  a proof,  see  [18]  and  [19].) 

We  note  that,  for  ease  of  exposition,  we  have  worked  over  the 
quotient  field  K of  R.  However,  if  the  data  were  defined  over  R, 
then  the  realization  problem  is  essentially  a question  of  rationality 
which  accounts  for  the  assumption  of  integral  closure.  Thus,  the 
results  presented  here  hold  with  data  in  R (see,  (43],  [44],  (47]). 


Classical  Groups  of  Symmetry  and  Moduli  for  Linear  Systems 
There  are  three  main  topics  covered  here.  The  first  is  the 
construction,  in  external  and  in  internal  terms,  of  the  moduli  space 
for  linear  systems  and  a proof  of  the  theorem  on  "realization  with 
parameters",  which  is  based  on  Zariski's  Main  Theorem.  The  second 
is  a brief  survey  of  biregular  and  topological  invariants  of  the 
moduli  space,  with  emphasis  on  the  impact  these  invariants  have  on 
system  theory,  as  in  Section  2.  We  close  by  considering  moduli 
for  systems  having  richer  external  symmetries,  such  as  Hamiltonian 
systems  or  systems  describing  symmetric  electrical  networks.  Here 
we  see  a more  arithmetic  side  come  to  the  foreground,  and  in  fact 
we  sketch  the  equivalence  between  an  old  problem  in  network  synthesis 
and  what  might  be  called  the  "Serre  conjecture  for  quadratic  modules". 
In  particular,  recent  results  on  quadratic  modules  ([1])  turnout 
to  be  sufficient  for  the  classical  formulation  of  this  problem. 

Now,  over  JR  or  C,  it  is  natural  to  view  the  set  of  transfer 
functions  with  m inputs,  p outputs,  and  fixed  McMillan  degree  n 
as  a (complete)  metric  space.  Indeed,  this  topology  ([13])  is 
simply  the  compact-open  topology  on  transfer  functions,  regarded 
as  maps 

Tct:  JP1  -*■  Grass(m,m+p)  , 

as  in  (3.4).  Although  such  a description  is  in  fact  useful  in 
certain  contexts,  it  is  not  entirely  clear  that  this  space  is  an 
algebraic  variety,  nor  is  it  even  clear  what  its  dimension  should  be. 
For  this  reason,  the  first  problem  we  consider  is  to  construct,  as 
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explicitly  as  possible,  a moduli  space  for  the  action  of  GL(n,k) 

on  Sn  . In  fact,  since  it  is  rather  important  as  a classifying 
p ,m 

space  for  families  of  linear  systems,  it  is  rather  crucial  that  we 

construct  a fine  moduli  space.  This  is  done  in  several  stages,  it 

turns  out  (by  universality)  to  be  sufficient  to  construct  a geometric 

quotient  for  the  action  of  GL(n,k)  on  controllable  pairs  (A,B) 

2 

in  A1  . This  is  done  within  the  context  of  geometric  invariant 
theory,  the  rather  remarkable  result  being  that  the  properly  stable 
points  in  lAn  +nm,  relative  to  the  standard  character  linearization, 
is  precisely  the  space  of  controllable  pairs! . We  start  birationally . 

First  of  all,  the  corresponding  space  H of  Hankel's  is  clearly 
a quasi-affine  variety,  contained  in  Ifi^ , for  N = mp(n+l).  This 
was  shown  to  be  non-singular  by  J.M.C.  Clark  ([16])  and,  in  addition, 
his  proof  shows  that  IH  is  a special  variety,  in  the  sense  of 
Chevalley.  In  particular,  H is  rational  and  one  can  easily  write 
down  generators  for  the  function  field  . On  the  other  hand,  if 
R = k[A,B,C] , K the  field  of  fractions  of  R,  then  Falb  has  shown 
in  [18] : 


Proposition  4.1, 


— 1H  " * 


GL(n,k) 


This  follows  from  the  statement  that  the  ring  of  invariant 
regular  functions  defined  on  the  Zariski  open  subspace  of  semi-simple 
A is  generated  by  the  entries  of  the  Hankel.  This  is  proved, 


The  birational  situation  using  these  ideas  was  sketched  in  [18]  and 
the  entire  theory  fully  developed  in  [10] . 


for  B and  C vectors,  by  Gurevitch  - as  an  illustration  of  the 

symbolic  method!  It  is  proved  in  general  in  [18] . 

The  biregular  situation  involves  a bit  more  care  - this  is 

due  to  several  people  (e.g.  [14],  [26])  and  our  treatment  follows 

([14]),  relying  on  [40].  It  can  be  shown  that,  if  Zn  is  the 
- 2+ 

Zariski  open  subspace  of  /An  +nm  of  controllable  pairs  and  if  Rq 

is  the  ring  of  regular  functions  on  Zq,  then  in  contrast  to 

Proposition  5.1,  RpL^n’— does  not  separate  orbits  in  Zq,  unless 

m = 1.  In  fact,  if  k is  algebraically  closed, 

rGL (n , k)  = k[Ci) , (4.2) 

where  (c^)  denote  the  characteristic  coefficients  of  A.  In 
particular,  although  Theorem  3.9  implies  that  the  orbits  in  Z^ 
are  closed,  an  affine  quotient  does  not  exist,  unless  m = 1. 

However,  we  do  have: 

Theorem  4.3.  ZQ  -*■  ZQ  is  a universal  geometric  quotient,  with  ZQ 
a quasiproj ective , non- singular , irreducible  variety  of 
dimension  nm. 

Sketch  of  Proof.  Following  [40],  since  (4.2)  implies  that  the 
set  of  properly  stable  points  for  the  trivial  character 
(Pic(ZQ)  = (0)!)  is  empty,  we  compute  the  properly  stable  points 
for  the  relative  invariants  of  weight  1.  Here  we  find  a surprisingly 
pleasant  result: 
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Proposition  4.4.  The  set  of  properly  stable  points,  relative  to 


det(-)>  in  /An  +nm  is  exactly 


This  is  very  much  related  to  the  fact  that  Eq  is  the  principal 
orbit  type  for  the  action.  In  the  scalar  input  case,  this  can  be 
proved  by  linear  algebra,  one  direction  being  already  implied  by 
the  proof  of  the  state-space  isomorphism  theorem.  For  the  other 
half,  to  say  (A,B)  is  in  the  principal  orbit  type  is  to  say,  by 
differentiating  the  action, 


n = dim  0^  ^ = dim  0^  + dim  Z(A)B, 


(4.5) 


where  Z(A)  is  the  centralizer  of  A and  0 denotes  the  orbit 
of  • . Now,  (4.5)  can  be  refined  since 


dim  0.  < n - n and  dim  Z(A)B  < n 


(4.6) 


holds  for  any  A.  Thus,  dim  0^  = nz  - n is  maximal  and  A there- 
fore possesses  a cyclic  vector.  By  Frobenius's  Theorem,  Z(A)  = k [ A ] 
and  thus 


dim  Z (A) B = dim  k [A] B = rank(B,  AB , . . . , An" 1B)  = n, 


which  was  to  be  shown.  The  multivariable  case  also  follows  from 
linear  algebra  (see  D.4]  ) . 

Assuming  Proposition  4.4,  Theorem  4.3  follows  from  Mumford's 
Theorem  (ftOI,  Theorem  1.19). 
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By  universality,  a geometric  quotient  for  the  twisted  action 

of  Gi.(n,k)  on  Ep  = {(A,B,C):  (A,B)  controllable}  exists  and 

shares  many  of  the  same  properties.  In  fact,  Ep  -*■  Eg  is  an 

algebraic  vector  bundle  with  base  Eg  and  fiber  kn^ . Since  the 

quotient,  Ep  -*■  Ep  is  geometric,  we  also  obtain  a quotient  Ep  m 

for  the  Zariski  open  subspace  Sn  „ C £ . 

* P,m  P 


Theorem  4.7.  A fine  moduli  space  for  state-space  equivalence  in 

Sn  exists.  Moreover,  En  is  a rational,  smooth,  irreducible 
p,m  p,m 

quasiproj ective  variety  of  dimension  n(m+p) . 


The  construction  of  a universal  family  of  systems  on  Eg 
was  carried  out  in  [26]  and  extended  to  Ep  m in  [10].  Thus,  any 
map  f:  X ■+  Ep  m induces  a linear  system  defined  over  the 
coordinate  ring  R of  X.  This,  of  course,  extends  to  the  study 
of  differentiable  or  continuous  families  as  well. 


Theorem  4.8.  If  H(x)  is  a regular  Hankel  matrix,  of  locally 
constant  rank  on  X,  then  H(x)  is  realizable  over  R,  with 
finitely  generated,  projective  state  module. 


Proof.  After  all  of  the  above,  including  the  existence  of  a 
universal  family,  i.e.  a vector  bundle  Q - the  state  bundle,  an 
endomorphism  srf  of  Q,  and  sections  and  co-sections  j , 

resp.,  this  amounts  to  proving  that  realization  is,  in  fact, 
algebraic  in  its  parameters  (see  [25],  [14]).  That  is,  if  n is 
the  natural  map 
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n:  E®  m + H (4.9) 

p,m 

[A,B,C]  - (CAi+^  _1B) 

we  claim  n is  biregular.  By  Theorem  3.9,  n is  bijective  and 
by  Proposition  4.1  n is  birational.  By  Clark's  Theorem,  M is 
normal  and  thus  n is  biregular  by  Zariski's  Main  Theorem.  Finally, 
the  assumption  on  H(x)  is  simply  that  H define  a map  to  IH. 

Example  4.10.  Given  (h^)  e ^ (^) ^n+1^  of  rank  n,  we  seek  a 
realization  (A(z) ,B(z) ,C(z) ) summable  in  z e 2.  If  we  ask  that 
the  dual  system  be  minimal  too,  the  following  sufficient  condition 
(see  [10],  [46])  becomes  necessary. 

(T)  the  closed  linear  span  of  det(huj)  and  its  translates 
is  i,1  (2)  . 

Indeed,  Fourier  transforming  and  using  the  Tauberian  Theorem, 
we  obtain  a continuous  map 

Ch^)  : S1  - H (4.11) 

and,  thus,  a linear  system  (A,B,C)  defined  over  S,1^).  since  n 1 
is  analytic.  Now,  the  state  module  is  trivial  over  CCS1),  for 
topological  reasons,  and  the  Docquier-Grauert  Theorem  asserts  that 
the  corresponding  bundle  can  even  be  trivialized  over  S,1^)!  In 
particular,  we  obtain  a realization,  in  matrix  form,  by  inverting 
the  Fourier  transform.  Recent  work  by  M.  Hazewinkel  on  constructing 
partial  compactif ications  for  m,  suggests  that  we  may  be  able  to 
drop  the  hypothesis  on  the  Hankel . 
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r 

I t 

i ! 
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We  stress  the  fact  that  the  computation  made  above  actually 

depends  on  the  topology  of  S^;  for  example,  it  can  be  shown  that 

there  exists  IP^  C En  , of  codimension  nfp+1) , over  which  the 

p,m 

universal  family  is  non-trivial.  This  is,  of  course,  very  much 
related  to  the  existence  of  algebraic  (or  continuous)  canonical 
forms  - a problem  which  has  received  a great  deal  of  attention. 

Now,  Mumford's  Theorem  yields  an  actual  imbedding  of  Eg  in 
¥ ; i.e.  a very  ample  line  bundle  L on  Eg.  L also  arises  in 
another  way,  viz.  L = AnE^,  as  one  can  see  by  checking  cocycles 
(c.f.  [14],  Section  4),  and  thus  L is  associated  to  the  principal 
GL(n,k) -bundle,  Eg  -►  Eg.  Furthermore,  L cannot  be  trivial  if 
m > 1,  by  (4.2).  If  m = 1,  Eg  = /An  by  the  existence  of  the 
rational  canonical  form.  The  more  precise  results  for  m > 1 are: 


Theorem  4.12.  If  k is  algebraically  closed,  Pic(Eg)  = (L)  - 2 
and  there  exist  no  algebraic  canonical  forms.  If  k = C,  then  the 
Betti  numbers  do  not  vanish,  for  i = 0,...,nm  - n and 

continuous  canonical  forms  do  not  exist. 

Sketch  of  Proof.  Since  Eg  -*•  Eg  is  an  algebraic  principal 
GL(n,k) -bundle,  one  can  conjugate  Pic(EQ)  by  descent  ([40],  p.  32): 

picGL(n,k)(£o)  s pic(Eo). 

Since  Pic(Eg)  = (0),  we  have  a surjection  (by  the  rigidity  lemma), 


Lt 


GL (n,k) 


PicGL (n, k) 


t 


[ 

( 

i. 
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which  is  injective  by  (4.2).  In  terms  of  group  cohomology,  this 
is  just  the  map, 

GL(n^k)  - H1(GL(n,k)  ,k*)  - H^GLCn.k)  ,0*  ) = Pic61^11*^  (In), 

- - - Lo 

I 

induced  by  the  inclusion  of  GL  (n, k) -modules , k . If  m = 1, 

the  unit  det (B , AB , . . . , An _1B)  is  a coboundary  of  the  character 

0 ,0  * 

det(*)  , regarded  as  a cocycle  with  value  in  47  , 

0 

As  for  the  topological  obstructions,  notice  that  (4.2)  asserts 
that  IAn  is  the  categorical  quotient  for  the  action,  in  the 
category  of  affine  varieties.  In  particular,  consider  the  natural 
foliation 


In  the  codimension  1 case,  n = 1 and  thus  Zq  = * IP™  1 while 

Zq  -*•  Zg  is  the  pullback  of  the  universal  fibration.  This  is  quite 
general,  on  a Zariski  open  subspace  of  /An,  the  leaves  are  co- 
dimension n projective  spaces,  IPN!  Over  (E,  these  can  be  imbedded 
non-trivially , giving  a proof  of  the  second  statement.  Q.E.D. 

Over  IR,  M.  Hazewinkel  ([24))  has  shown  that  there  exist  no 
continuous  canonical  forms,  by  a clever  imbedding  of  S1  into  ZQ. 
It  is  possible  to  prove  non- vanishing  of  higher  homotopy  groups 
in  the  proper  ranges  over  IR  and  <C,  by  making  some  crude 
computations  along  the  lines  in  [11] , and  applying  Bott  Periodicity. 
These  computations  also  show  that  the  natural  foliation  on  Zq  is 
not,  in  general,  a product  of  /\n  with  a typical  leaf.  Of  course, 
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the  situation  remains  the  same  for  I , but  En  is  a bit  more 

p P,m 

elusive.  However,  the  statements  about  canonical  forms  still  hold, 
answering  in  the  negative  - by  topological  and  geometric  means  - a 
question  which  many  engineers  had  assumed  to  be  evident.  This  global 
form  of  the  question  was  originally  raised  by  R.  E.  Kalman  [31]. 

For  the  remainder  of  the  section,  we  consider  the  real  field 
and  discuss  an  application,  following  [13],  of  a theorem  of 
Harder  to  a problem  in  network  synthesis,  posed  by  Youla  in  [53]. 

It  should  be  stressed  that,  even  a decade  ago,  this  problem  was 
known  to  be  algebro-geometric  in  nature  and,  in  fact,  the  partial 
results  obtained  by  Koga  ([36])  were  based  on  his  use  of  the 
Riemann-Roch  theorem  for  curves  to  analyze  coprime  factorizations  of 
transfer  functions  in  two  variables. 

Now,  a result  of  Glover  ([21])  asserts  that  if  m > 1 or  p > 1, 
the  space  H is  connected  in  contrast  to  Theorem  2.7.  This  is  due 
largely  to  the  fact  that,  even  in  case  m = p,  the  Hankel  is  not 
necessarily  symmetric  and  one  has  more  room  to  move  in.  However, 
symmetric  Hankels  or,  equivalently,  symmetric  transfer  functions 
are  important  in  network  theory  (see  [6]  , [37],  [54]).  This 
symmetry  is  displayed  in  state  space-form  as: 
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([4])  as  the  signature  of  the  indefinite  metric  in  the  Brayton-Moser 
equations.  In  fact,  the  state-space  isomorphism  theorem  applies  here 
to  give  the  group  of  symmetrices  ^(p,q)  C GL(n,]R)  (see  [54]). 

Moduli  for  these  linear  dynamical  systems  have  been  studied  in  [13] 
where,  in  particular,  Theorem  2.7  is  shown  to  hold.  The  invariant 
playing  the  role  of  the  Cauchy  index  is  the  Maslov  index  where  the 
real  Grassmann  in  (3.4)  is  replaced  by  a Lagrangian  Grassmann. 

Now,  the  problem  alluded  to  above  is  just  a question  of 
realization  with  parameters,  subject  to  the  symmetry  constraints. 

That  is,  given  a symmetric  transfer  function  involving  delays  - 
which  is  simply  a transfer  function  over  R[x^,...,x^]  (.[34],  [35])  - 
find  a symmetric  realization.  By  the  validity  of  the  Serre 
conjecture,  the  state  module  can  be  assumed  free  and  the  corre- 
sponding question  is  equivalent  ([13])  to: 


given  a symmetric  Hankel  with  constant  rank,  defined  over 
R =lR[x],  find  a symmetric  realization  (4.13)  over  R. 

Since,  over  1R,  this  depends  only  on  Sylvester's  Theorem,  we 
can  invoke  the  following  to  finish  off  the  problem. 

Theorem  (Harder  [1]).  Any  non-degenerate  symmetric  bilinear  form 
defined  over  k[x]  is  equivalent  to  a form  defined  over  k. 

We  close  by  remarking  that  the  general  theory  of  quadratic 
modules  arises  in  several  instances,  especially  in  connection  with 

families  of  Hamiltonian  systems.  Here,  again,  over  1R  we  have  a 
Hamiltonian  realization,  with  symmetry  group  Sp(.n^R)  and  thus  the 
moduli  problem  is  already  "in  canonical  form". 

lI  


A 


5 . Stabilization  and  Feedback  Groups 

To  be  sure,  one  of  the  most  important  results  in  deterministic 
control  theory  is  the  pole-placement  theorem,  Theorem  5.2,  which  is 
concerned  with  stabilization  of  linear  systems  by  means  of  "state 
feedback"  - a topic  very  much  related  to  the  material  developed  in 
Section  2.  It  is  fair  to  add,  however,  that  a much  deeper  under- 
standing of  the  whole  feedback  question  comes  only  with  the  interpre- 
tation of  feedback  substitutions  as  part  of  a group  action  on  IA^ , 

2 

N = n + nm  + np.  Explicitly,  for  (n,m)  fixed  we  define 
9 = ^(n,m)  as  the  semi-direct  product  of  3 groups,  generated  by 
the  actions  on  triples  (A, B,C): 


(i)  GL(n,k)  acting  via  change  of  basis  in  the  state  space, 

(ii)  GL(m,k)  acting  via  change  of  basis  in  the  input  space, 

(iii)  Hom(kn,km)  acting  via  (A,B)  h-  (A+BF,B). 

This  action  is  more  complicated  than  the  action  in  (i),  for  both 


mathematical  and  control-theoretic  reasons.  That  is,  9 is  a non- 
reductive  extension  of  GL(n,k)  and,  although  9 preserves  con- 
controllability,  and  hence  acts  on  9 does  not  leave  the  space 

of  linear  systems,  Sp  m,  invariant.  Thus,  we  begin  this  section  by 
studying  the  action 


9 * I0  -*  IQ  (5.1) 

and  extending  to  the  case  p > 0.  This  extension  is  possible,  in 
the  sense  of  constructing  moduli,  since  each  of  the  isotropy  sub- 
groups arising  in  (5.1)  is  the  semidirect  product  of  a reductive 
group  with  the  unipotent  radical.  We  reproduce  the  theorem  of 
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Hermann  and  Martin  ([28])  which  states  that  the  classical  invariants 
for  (5.1)  - which  are  arithmetic , (being  the  Kronecker  set  introduced 
in  Section  3)  - arise  as  the  Birkhof f -Grothendieck  invariants  of  a 
bundle  on  P*.  We  extend  these  considerations  to  families  of 
controllable  pairs  and,  following  Harder,  Narasimhan  and  Shatz  ([45]), 
introduce  what  might  be  called  the  "ubiquitous  ordering"  on  the  set 
of  partitions  of  an  integer  n.  Indeed,  this  ordering  is  well-known 
in  system  theory,  having  been  introduced  by  Rosenbrock  ([42])  in  his 
study  of  the  generalized  pole-placement  theorem,  and  we  briefly 
indicate  (based  on  joint  work  with  C.  Martin)  how  this  circle  of 
ideas  yields  some  powerful  new  tools  in  the  analysis  of  state 
feedback. 

First  of  all,  recall  that  the  poles  of  a transfer  function  T(s) 
are  precisely  the  roots  of  xA  - the  characteristic  polynomial  of  A, 
where  (A,B,C)  is  a linear  system  realizing  T(s).  Now,  if 
F e Hom(Rn,Rm),  (R  a commutative  ring  with  identity)  we  may  regard 
the  action  (iii)  of  F as  feeding  back  the  state  of  the  system  as 
a component  of  the  input.  In  this  setting,  we  ask  (over  ]R)  whether 
there  exists  an  F such  that  XAj.t>c  is  a Hurwitz  polynomial. 

More  generally,  over  an  arbitrary  field  k we  consider  ([51]): 

Theorem  5.2.  Any  monic  p(s)  can  be  expressed  as  xA+Bp  if>  and 
only  if,  (A,B)  is  controllable. 

In  the  scalar-input  case,  this  result  was  classically  derived 
from  the  existence  of  the  rational  canonical  form,  and  thus  holds  for 
any  ring  R.  However,  for  families  of  systems  with  m > 1, 
technicalities  arising  from  the  non-existence  of  canonical  forms, 
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Theorem  4.12,  are  non-trivial.  In  fact,  a counterexample  for  the 
general  situation  has  recently  been  announced  by  Bumby  and  Sontag 
([9]).  In  general,  our  knowledge  for  arbitrary  coefficient  domains 
is  very  scant  and  is  complete  only  under  strong  assumptions  on  the 
ring  and  the  type  of  monic  polynomial  considered  ([39])  or  on  the 
type  of  system  ([12]). 

Nevertheless,  we  can  still  glean  some  insight  from  the  scalar 
input  case.  In  terms  of  the  action  (5.1),  which  is  transitive  in 

this  case,  we  fix  (A,B)  and  consider 


a (A , B)  * * * R ’ 


(5.3) 


“(A,B) (g’s’F)  = Cc± (g (A+BF) g )) 


(where  c^(-)  is  the  i characteristic  coefficient).  Following 


E.  Kamen , we  note  that  the  image  of  a 


(A,B) 


is  the  translation, 


by  the  vector  xA»  of  a linear  subspace,  so  that  if  we  reduce 
modulo  (xA),  and  supress  g and  s,  the  corresponding  map 


a ( A , b ) " Hom(Rn,R)  Rn 


(5.4) 


is  linear!  If  ^(A,B)  _ [u>AB,...,An  ^B]  is  the  controllability 
operator,  then  the  relationship  between  if  and  a,  i.e.  Theorem  5.2, 
is  brought  out  by  the  beautiful  identity  due  to  B.F.  Wyman  ([52]): 


coker  a = coker (if  ) 


(5.5) 
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If  the  pair  (A,B)  is  controllable  over  the  fraction  field  of  R, 
the  right  hand  side  can  be  computed  in  terms  of  % , i.e.: 

Ext  * (coker  5f,R)  = coker  a,  Ext  *(coker  a,R)  = coker  Sf.  (5.6) 

In  particular,  if  R is  a field  and  if  A is  semi-simple,  (5.5) 
contains  the  folklore  result  that  one  cannot  change  the  uncontrollable 
modes  by  state  feedback.  That  is,  if  A is  diagonal,  then  any 
eigenvalue  A^  for  which  b^  is  zero  must  persist  in 

In  the  case,  m > 1,  if  k is  algebraically  closed  and  of 
characteristic  0,  we  can  replace  a by  its  Jacobian  and  thus  obtain 
the  corresponding  generic  result  from  the  dominant  morphism  theorem. 
The  computation  of  the  Jacobian,  over  (E,  is  due  to  Hermann  and  Martin 
([27])  and  relies  on  a differential  analogue  of  Proposition  4.4. 

Thus  far,  the  only  invariants  which  we  have  attached  to  a 
linear  system  have  been  topological,  viz.  the  degree  of  the  associated 
transfer  function, 

T:  IP^  -+  Grass^(m,m+p)  . (5.7) 

For  k = 1R,  this  defines  the  Cauchy  index  if  m = p = 1 or  else 

the  McMillan  degree  mod  (2).  For  k = C,  this  defines  the  McMillan 

degree.  However,  T is  an  algebraic  map,  and  as  such  determines 

* 

finer  invariants,  namely:  the  isomorphism  class  of  T U,  where  U 
is  the  (algebro-geometric)  universal  vector  bundle  on  Grass^(m,m+p) . 

By  the  well-known  result  of  Birkhoff — Grothendieck  (and  probably 

* 

others!),  T U is  determined  by  a partition  of  deg^T  into 
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I 

non-negative  integers  (n^)  , ni  . n2  - * • • - nm  > 0 , if  k is 

algebraically  closed.  Clf  k = ]R,  a theorem  of  Serre's  still  gives  a 

* 

decomposition  of  T U into  a sum  of  line  bundles,  but  these  in  turn 
are  only  determined  by  integers  mod(2).)  We  owe  to  Hermann  and 
Martin  the  identification  of  these  geometric  invariants  ([28]): 

* 

Theorem  5.8.  The  Birkhoff -Grothendieck  invariants  of  T U are  the 
Kronecker  set  of  any  minimal  realization. 

Sketch  of  Proof.  Since  one  set  of  invariants  is  defined  in  the 
frequency  domain  and  the  other  via  state-space  techniques,  the  key  step 
will  involve  the  Laplace  transform.  Note  that,  for  constant  coefficient 
differential  operators,  the  Laplace  transform  has  its  range  in  a field 
of  transcendence  degree  1 and  genus  0.  This  theorem  very  clearly 
indicates  the  relevance  of  the  genus  in  such  questions. 

The  key  step  is  due  to  Rosenbrock:  over  <E,  consider  the  Laplace 
transforms  of  the  solutions  to  the  initial  value  problem, 

*(t)  = Ax(t)  + Bu(t)  , x(0)  = 0. 

These  are  elements  of  the  <C[s]  -module,  defined  as  the  kernel  of 
the  following  pencil  of  matrices, 

[A-sI  ,B]  : (<En  © <Cm)  ® <C[s]  <Cn  ® C[s]  . (5.9) 

The  theory  of  such  pencils,  particularly  the  problem  of  equivalence, 
is  due  to  Weierstrass  and  Kronecker  (see  [20])  and  applies 
remarkably  well  to  the  situation  at  hand.  Explicitly,  controllability 

d 
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ba* 

of  (A,B)  is  equivalent  to  right  surj ectivity  of  (5.9)  and  the  notion 
of  strict  equivalence  of  pencils  is  precisely  that  of  feedback 
equivalence  of  (A,B)  with  (A',B').  Indeed,  strict  equivalence 
is  an  equivalence  of  k[s] -module  maps,  with  the  isomorphisms 
independent  of  s.  By  writing  the  definition  of  strict  equivalence 
and  comparing  degrees,  we  find  that  the  group  of  strict  equivalence 
of  pencils  of  the  form  (5.9)  is  a matrix  representation  of  &\ 

Since  independence  of  s connotes  proj ectivity,  we  follow  another 

clue  from  Kronecker:  we  projectivize  the  pencil  (5.9).  In  this 

l I 

setting,  right  surjectivity  guarantees  that  we  obtain  a homogeneous 

2 1 
bundle  on  A {0},  and  hence  an  m-bundle  on  1P^.  According  to  the 

Kronecker  recipe,  given  at  the  end  of  Section  3,  this  bundle 

m m 

decomposes  into  a sum  © <^(  k. ) , l <.  = n.  That  this  bundle, 

i=l  1 i=l  1 
1 * 

ker [tA-sI , tB] , on  IP^  is  isomorphic  to  T U,  follows  at  once  from 
observability:  we  simply  map  [x,u]  to  [u,Cx],  sending  the  kernel 

above  to  the  graph  of  T.  Observability  is  then  just  the  statement 
that  this  map  is  injective. 

Thus,  Kronecker's  classification  also  shows  more,  as  noted 
in  [32],  for  controllable  pairs: 

Theorem  5.10.  (A,B)  = (A',B')  mod  & if,  and  only  if,  the 

Kronecker  sets  are  the  same. 

In  fact,  Kronecker's  proof  (see  [20])  is  to  construct  a 
canonical  form  for  the  pencil,  given  the  data  ( k4).  This 
canonical  form  coincides  with  the  form  discovered  by  P.  Brunovsky 
([7])  by  other  methods.  We  remark,  however,  that  as  an  immediate 
application  of  the  construction  of  this  form,  pole-placement 

* — k 
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becomes  trivial.  Indeed,  the  invariants  n^  > n^  > •••  > > 0 are 

much  more  refined,  they  give  a "generalized  pole-placement  theorem". 


Theorem  5 . 11 . (Rosenbrock  [42]).  The  invariant  factors  of  A + BF 
can  be  made  arbitrary,  subject  only  to  the  constraints. 


A/ 

deg  > n1,  deg  <P  2 > n1  + n2,...,deg 


Another  way  of  interpreting  T U as  a complete  feedback 
invariant  for  the  action  (5.1)  is  to  count  sections.  Thus,  if 
h^(V)  is  the  dimension  of  H^QP^;V),  we  have 


h° (T*U)  = h°(®<?(n,))  = ? (n.+l)  = n + m, 

1 i = l 1 


(5.12) 


On  the  other  hand,  evaluating  T U at  “ gives,  canonically,  the 

& A 

input  space.  Thus,  if  <P : T^U  ■+  T2U  is  an  isomorphism,  <P  in- 

0 1 * 0 1 * 

duces  an  isomorphism  <P*:  H QP  ;T^U)  ^ H (IP  ;T2U),  commuting  with 
evaluations  at  00 ! Thus,  9*  is  a triangular  element  of  GL(n+m,k), 
corresponding  to  the  matrix  representation  of  <P*  obtained  above. 

This  point  of  view  has  several  applications.  For  instance,  if 
(A,B)  e Eg,  then  ^{k  B)  " t^ie  isotroPy  subgroup  for  the  action 
(5.1)  - is  just  Aut(V),  where  V is  the  bundle  defined  by  the 
kernel  of  [tA-sI,tB].  In  particular, 


dim  9/ 9 * dim  9 - h°(V  <g>  V*) 


dim  9 - ( l (n. -n.+l)) , 
ni-nj  3 


(5.13) 


k 
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a formula  originally  derived  by  R.  Brockett  in  [5],  by  a more 
explicit  representation  o f ^(A.B)  together  with  a heuristic 
argument  using  Young  diagrams  corresponding  to  the  partitions. 

Based  on  this  description,  it  can  also  be  shown  ([5])  that 
connectivity  Cover  ]R)  of  the  orbits  is  determined  by  the  Kronecker 
indices  mod(2). 

Now,  g-j  has  also  been  written  down  explicitly  by  Falb 

and  Wolovich  in  [19],  in  order  to  study  the  action  of  & (and  an 
extension  of  on  Ep.  Then,  one  studies,  in  addition  to  the 

Kronecker  set  of  (A,B),  the  action  of  on  An^.  In  [19], 

it  is  shown  that  a moduli  space  for  this  action  exists,  using  the 
fact  that  B)  *s  a^waYs  t*ie  semi-direct  product  of  a re- 

ductive group  with  its  unipotent  radical. 

For  the  remainder  of  the  section  we  consider  only  & acting 

on  controllable  pairs  CA,B)  and,  as  in  Section  4,  we  begin  to 

introduce  parameters.  Thus,  we  suppose  CA,B)  is  an  algebraic 

family,  defined  on  an  irreducible  affine  variety  X.  If  (A,B) 

arises  from  a delay-differential  system,  it  can  be  shown  that  we 
N 

can  take  X = Ik  ([10])  and  that  equivalence  modulo  ^(R)  is  not 
just  a formal  equivalence  ([35],  [10]).  In  this  case,  one  may 
assign  geometric  data  to  (A,B),  viz. 


V (A , B ) - X * Pk’ 

N 

as  before,  by  homogenizing  the  Rosenbrock  pencil.  If  X * A , it 
may  be  shown  ([10]),  based  on  the  alternative  view  presented  above, 


that 


(A,B) 


is  actually  a complete  invariant  for  (A,B)  under  ^(R) . 
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This  situation  is  an  improvement,  as  quite  a bit  more  is  known  about 
deformations  of  bundles  on  IP'*'  than  about  pencils  over  R.  For 
example,  if  the  deformation  is  point-wise  trivial,  i.e.  if  the 
Birkhof f -Grothendiech  invariants  are  constant  in  X,  then  a 
theorem  of  C.  Hanna  ([23])  asserts: 

V » Epjcv.)  « pJcV  C5'14) 

where  the  are  bundles  on  X,  and  the  VT  are  bundles  on  P"*- . 

N 

Thus,  if  X = A , such  a deformation  is  globally  trivial,  since 
is  necessarily  trivial.  Applying  this  to  ^ , we  find: 

Proposition  5.15.  If  the  pointwise  Kronecker  invariants  are  constant, 
CA,B)  is  feedback  equivalent  to  a system  defined  over  k. 

This  has  the  immediate  corollary  that  such  a system  is  coefficient 
assignable,  i.e.  Theorem  5.2  holds.  This  latter  fact  generalizes 
to  arbitrary  R,  although  the  Proposition  fails  to  hold  in  general 
([12]).  It  is  not  hard  to  show  that,  in  this  case,  the  pointwise 
Kronecker  invariants  are  equal  to  the  "global"  Kronecker  invariants, 
i.e.  the  Kronecker  invariants  computed  over  the  fraction  field  of  R. 
Indeed,  in  general,  there  is  a Zariski  open  set  on  which  the  pointwise 
invariants  are  constant  and  equal  to  the  global  Kronecker  in- 
variants. What  is  more,  the  exceptional  values  must  bear  some  re- 
in 

lation  to  the  global  invariants.  Explicitly,  if  V.  = © ^(n.), 

1 * i 1 

m f 1=1 

m , 

V'  = © <?(n!)  are  bundles  on  P , then  V > V*  in  the  Harder- 

i-1  1 

Narasimhan  ordering  just  in  case 
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| 

| 


1 i 

Z n.  > Z n!  , for  £ = 1 , . . . ,max(m,m' ) . (5.15) 

i=l  1 ' i=l  1 

Theorem  5.16  (Shatz,  [45])  The  Harder-Narasimhan  ordering  is  upper 
semi-continuous  in  algebraic  families  of  vector  bundles  on  p\ 

This  theorem  has  several  applications  in  problems  of  computing  feed- 
back invariants.  Now,  any  algebraic  family  of  controllable  systems 
induces  a map, 


and  Theorem  5.16  thus  gives  a relation  between  the  quotient  topology 
on  the  finite  set  £ & and  the  Harder-Narasimhan  ordering  or, 
equivalently,  the  Rosenbrock  ordering  of  Theorem  5.11.  We  call  the 
ordering  induced  by  specialization,  the  geometric  ordering,  while  we 
refer  to  the  "ubiquitous  ordering"  as  the  natural  ordering.  The 
natural  ordering  is  also  present  in  combinatorics  and  in  number 
theory,  where  it  is  studied  as  an  ordering  on  - the  set  of  all 

partitions  of  the  integer  n.  Now,  it  is  an  unpublished  piece  of 
folklore,  not  too  hard  to  prove,  that  a stronger  version  of 
Theorem  5.16  holds  for  algebraic  families  constructed  via  system 
theory,  viz.  the  geometric  ordering  is  dual  to  the  natural  ordering. 
In,  this  context,  the  following  theorem  ([8])  is  rather  remarkable: 

Proposition  5.17.  </  is  a lattice  under  the  natural  ordering. 

We  close  by  indicating  extensions  of  this  situation  which  are 
motivated  by  control  theory  but  seem  likely  to  contribute  some  new 
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examples  to  geometry.  First,  it  is  of  much,  more  practical  significance 
to  study  the  action  of  the  output  feedback  group;  i.e.  Hom(kn,km) 
is  replaced  by  HomC.k^,km)  in  (iii)  . In  this  setting,  the  proper 
generalization  of  Theorem  5.2  is  not  even  conjectured  - although  a 
counter-example  to  the  dominant  morphism  approach  has  been  constructed 
by  J.C.  Williems  ([48]).  Second,  pencils  may  also  be  defined  by 
control  systems  governed  by  partial  differential  equations.  In  this 
setting,  the  generalization  of  the  Rosenbrock  pencil  leads  to  coherent 
sheaves  on  F but  whether  these  all  split,  as  in  Grothendieck ' s 
theorem,  is  unknown  at  present.  Finally,  there  are  several  computations 
which  suggest  that,  for  non-constant  coefficient  differential  equations, 
the  transfer  function  ought  to  be,  in  the  scalar  input-output  case, 
a meromorphic  function  on  a Riemann  surface  of  higher  genus.  Each 
of  these  questions  deserves  further  study. 
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